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A mesoscopic Tera-hertz pulse detector 
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We show that, under the passage of an electromagnetic terahertz pulse, an asymmetric double 
barrier device may act as an on/off current switch, depending on the bias. The time-dependent 
response of the device is discussed. 
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It is well known that the characteristic frequencies of 
electronic processees in mesoscopic systems are in the ter- 
ahertz range. This follows from transition energies that 
are typically in the meV scale. For example, instabilities 
in transport through asymmetric double barrier, systems 
(ADBS) may give rise to terahertz oscillations S'EI Other 
time dependent processes, such as the charging and dis- 
charging of the well in these systems, are also expected 
to take place in the range of picoseconds. This suggests 
that an ADBS may react as a fast switch to the passage 
of terahertz pulse, a possibility that we explore in this 
work. 

ADBS are characterized by a bistable region of the 
bias, produced by charge accumulation in the space be- 
tween the barriers (the well). El The collector barrier is 
made wider with the purpose of increasing the lifetime 
of the resonance in the well, thus enhancing the amount 
of charge that is retained when current goes throughD'El 
At a critical bias V c [ the current drops abruptly due to 
a sudden emptying of the well, driven by an instability 
that may be understood using a nonlinear model. u The 
dynamics is dominated by the accumulated charge which, 
in effect, lifts the bottom of the potential well thus re- 
taining the resonance condition beyond V c i for ballistic 
transmission of an incoming electron. There is a point at 
which this charge is so large that it becomes favourable 
to spill it out, the well is emptied and the resonance con- 
dition is lost, followed by a large current drop. Similarly, 
when the well is uncharged it will remain so, even as the 
bias is decreased to values smaller than V c i- At a sec- 
ond critical value V c -\ the resonance condition is fullfilcd 
again and current now flows. This completes the bistable 
cycle, a signature of which is the fact that V c i < V c -\ . 

In this work we propose that an ADBS device biased 
slightly below V c y will still undergo a transition, triggered 
by the passage of radiation in the terahertz region. The 
external field introduces an additional oscillating field 
that may in effect bring the bias to criticality. We also 
contend that a radiation field may trigger the onset of 
resonant transport if the system is biased slightly above 
Vc]- Optical radiation will have no effect in either case 
since the field then oscillates so fast that the electrons 
have no time to respond. Only at terahertz and lower 
frequencies would one expect the system to switch from a 



state of high current to one of low flow of electrons in . 
presence of an incident radiation pulse, or vice-versa.B 

Consider an ADBS under bias and in the presence of 
an electromagnetic field polarized along the z-axis, the 
growth direction. In order to study the time evolution of 
the device we adopt a first-neighbors tight-binding model 
for the hamiltonian. The radiation field enters as a space 
and time dependent voltage. To a good approximation 
the longitudinal degrees of freedom are decoupled from 
the transverse motion and may be treated independently. 
The probability amplitude b" for an electron in a time 
dependent state \a >, to be at pjLfne j along z, is deter- 
mined by the equation of motiorn 
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In this expression ej(t) includes the fixed band contour, 
the external radiation-induced voltage bEsin2ixvt and 
the applied dc bias, the latter represented by a term lin- 
ear in the spatial coordinate j. The sum over (3 covers 
all occupied electron states and v is the hopping matrix 
element between nearest neighbor planes. In writing Eq. 
(1) we have adopted a Hartree model for the electron- 
electron interaction, keeping just the intra-atomic terms 
as measured by the effective coupling constant UE As we 
will show in what follows, this nonlinear term is of key 
importance in the behavior of the system. 

The time dependent Eq.(l) is solved using a half- 
implicit numerical, method which is second-order accu- 
rate and unitaryntJ. Boundary conditions must be spec- 
ified at the left (z — —L) and right (z — L) edges of 
the structure. The approach taken here assumes that 
th||jy|Syc function al I hue / is »i\vu outside I he structure 
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Here k a and k' a = ^/2m*[e a — tj^fh are the wavenum- 
bers of the incoming and outgoing states, respectively, 
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with e a — — Av sin 2 (k a a/2) the energy of the incoming 
particle. To model the interaction with the particle reser- 
voir outside the structure the incident amplitude I is as- 
sumed to be a constant independent of the coordinates. 
The envelope function of the reflected and transmitted 
waves, Rj and Tj, are allowed to vary with j, however. 
Since far from the barriers these quantities are a weak 
function of the coordinate Zj we restrict ourselves to the 
linear corrections only. This approximation is appropri- 
ate provided the time step St does not exceed a certain 
limiting value. For the results presented here, a value of 
St = 3 x 10~ 17 s was found sufficient to eliminate spurious 
reflections at the boundary while maintaining numerical 
stability up to 40xl0~ 12 s. In our numerical procedure 
the coefficients obtained without electromagnetic field for 
a given dc bias are used as initial condition when the THz 
radiation is turned on. With &"(£) known, the time de- 
pendent current at site j is obtained numerically fromi 

Jj(t) = | J' Im{b* a (b« +1 - b«)}(kj ~ k 2 a )dk a , (4) 

where kj = ^/2m*e//?i, with e/ the Fermi energy. 

We next apply our model to an asymmetric 
GaAs/AlGaAs double barrier structure, with emitter and 
collector barrier thicknesess of 1.12 nm (2 sites) and 3.36 
nm, (6 sites) respectively, and a well thickness of 11.2 
nm (20 sites). The second barrier is made wider than 
the first in order to enhance the trapping of charge in 
the well. For this geometry the first resonance at zero 
bias occurs at 30 meV. The conduction band offset is 
set at 300 meV. The buffer layers are uniformly doped 
up to 3 nm from either barrier, so as to give a neutral- 
izing free carrier concentration of 2x 10 17 cm -3 at the 
contacts. In equilibrium, the Fermi level lies 19.2 meV 
above the asymptotic conduction band edge, so that the 
zero bias resonance lies well above the Fermi sea. The 
contribution to the potential due to the applied bias is 
taken into account throught a term linear in j, which is 
assumed to arise from fixed charges. The parameter val- 
ues in Eq. (1) are set at v — -2.16 eV and U = 100 meV. 
The latter was chosen phenomenologically so as to fit the 
experimental J-V characteristic for a GaAs devices.Q The 
sample has 400 sites and the normalization of the wave 
functions is chosen so that charge from the electrons fill- 
ing up to the Fermi eaergy exactly cancels the positive 
charge at the contacts.Q We solved Eq. (1) using the pro- 
cedure described above, for an energy mesh appropriate 
to compute the integral in Eq. (4) . Good convergence was 
found for a mesh of 100 points. 

Figure 1 shows the current-voltage characteristic in the 
absence (solid line) and presence (dashed line) of radia- 
tion of amplitude SE = 10 meV at v — 1 THz. In the 
latter case we exhibit an average of the current over time. 
Note that at the chosen values of parameters V c i = 0.320 
V and Vci = 0.282 V. It is clear from the figure that 
the radiation field narrows down the region of bistabil- 
ity, in agreement with previous results by Ifiarrea and 



PlatercO. This effect is to be expected since the time 
dependent field added to the bias brings the system pe- 
riodically to criticality when the applied dc bias has not 
reached this condition yet, thus triggering the charging 
or discharging of the well. 

In Fig. 2 we show the time evolution of the charge 
density at the center of the well for different frequencies 
of the radiation field at V = 0.310 V (empty circle in 
Fig. 1), a bias slightly below V c i- The condition for 
resonant tunneling is still met, conduction is allowed and 
the well is initially charged. A THz field of the same 
amplitude as for Fig. 1 is turned on at t = 0, and as the 
radiation passes through the system the well empties, 
doing so in a few picoseconds time. Two characteristic 
times are involved in the data: the external radiation 
period T r — \ jv and the time T w ~ 3 ps it would take 
for our well to empty if charge is initially in it. When 
T r 3> T w the radiation field essentially acts as an added 
dc bias and the well empties within the time T w , while 
in the other extreme T r <C T w the oscillation is so fast 
that the electrons cannot respond and the system remains 
charged and conducting. 

Figure 3 shows a situation in which the well is initially 
empty at a bias of V — 0.285 V, slightly above the crit- 
ical value V c \ (empty square in Fig. 1). It is physically 
reached by lowering the bias after it has gone beyond 
V c y. Once again the THz field is switched on at t = 0. 
We observe that in all cases exhibited the well begins to 
charge, and after a transient time the system enters full 
resonance and current flows. 

The above results assumed a radiation field of fixed 
amplitude SE = 10 meV. One may ask how close to the 
critical value V c ± must the system be biased in order to 
act as a switch for weaker radiation fields. This is shown 
in Fig. 4 for v = 0.33 THz and assuming the switching to 
take place at r = 17 ps time. The bias offset is defined as 
SV = Vci ~~ V ■ The region above the curve (labeled YES) 
is where the potential drop takes place within the time r, 
while the region below (labeled NO) is where the switch- 
ing does not take place in that time interval. Within the 
range of our calculations we found the shape of the curve 
in Fig. 4 to be generic, shifting upwards as the frecuency 
increases. Close to the origin the dependence is approx- 
imately linear and for the chosen frequency follows the 
relation SE ~ \{\ + SV). Using this expression we get 
that at a bias SV = 1 meV our device would switch under 
radiation of about 50 Watt/cm 2 and stronger. The sensi- 
tivity could be improved using a wider collector barrier, 
thus having a narrower resonance (longer lifetime T w ). 
Because of limitations due to numerical instabilities this 
ansatz would be best tested experimentally. 

In summary, we have shown that an asymmetric dou- 
ble barrier heterostructure may act as a switch triggered 
by the passage of electromagnetic radiation at frequencies 
in the terahertz region and below. The frequency thresh- 
old for this switching action depends on the barrier and 
enclosed well widths. Depending on the applied external 
bias, the passage of current is turned on or off by the ra- 
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diation pulse. Our results rely on the current drop as the 
resonance in the well falls below the emitter conduction 
band edge, a feature also present in symmetric double 
barrier heterostructures. In the latter case however, the 
drop is not an instability of the system and does not take 
place abruptly, a desirable feature for a switching device. 

Work supported by FONDECYT grants No. 1990425 
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FIG. 1. Current- Voltage characteristic for SE = (solid 
line) and SE =10 meV(dashed line) at v = 1 THz. 



FIG. 2. Time evolution of the charge density in the center 
of the well at V = 0.310V and frequencies 0.10 THz, 0.33 THz 
, 0.50 THz and 0.66THz. 



FIG. 3. Time evolution of the charge density in the center 
of the well at V = 0.285 V and frequencies 0.10 THz, 0.33 
THz , 0.50 THz and 0.66 THz. 

FIG. 4. Boundary separating the region at which switch- 
ing takes place within 17 ps of the arrival of an external pulse 
at frecuency v — 0.33 THz, from the region in which the well 
remains charged beyond that time interval. 
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